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Abstract 

Recent experimental studies have shown that well-annealed, unstrained Si(105) 
surfaces appear disordered and atomically rough when imaged using scanning tun- 
nelling microscopy (STM). We construct new models for the Si(105) surface that 
are based on single- and double-height steps separated by Si(001) terraces, and pro- 
pose that the observed surface disorder of Si(105) originates from the presence of 
several structural models with different atomic-scale features but similar energies. 
This degeneracy can be removed by applying compressive strains, a result that is 
consistent with recent observations of the structure of the Ge/Si(105) surface. 
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The self-organized growth of Ge/Si quantum dots has been investigated extensively 
for more than a decade, driven by their potential applications as optoelectronic devices and 
nanoscale memories. In the early stages of growth, the quantum dots that have pyramidal 
shapes bounded by {105} facets evolve from stepped mounds without encountering any ener- 
getic barriers for their nucleation ^ [2] . The absence of nucleation barriers has been explained 
by a competition between the strain-dependent, negative step formation energy and repulsive 
step-step interactions that have a weak dependence on strain [31 E|. The atomic configu- 
ration of the (105) facets has also been elucidated; a rebonded step model for the Ge(105) 
surface under mismatch strain was found to play a crucial role for the stability of this surface 

[3J[S1 © Ill- 
While the structure of Ge/Si(105) surface has been recently elucidated, the Si(105) 
surface shows intriguing features that are not well understood. Experimental work by Tomitori 
et al. 0, Fujikawa et al. (SJ and Zhao et al. [2] reveal that Si(105) is atomically rough even 
after careful annealing, and its STM image does not display two-dimensional periodicity. 
The analysis of Zhao and coworkers |2j suggests the presence of a structure for the Si(105) 
surface with large (001) facets and double-height steps. While a model for this surface with 
double-height steps was presented in a study of its stability has not been attempted. 

Motivated by the recent STM investigations Eh we search for reconstructions of 
Si(105) based on (001) terraces separated by single- and double-height steps oriented along 
(100) directions. In addition to the currently accepted models of Si(105) (shown, for example, 
in Ref. [31), we have found a few other possible structures based on double-height steps. We 
investigate the stability of these novel reconstructions, and propose that the roughness of the 
Si(105) surface is due to the coexistence of different structures with surface energies that are 
very close to each other. Furthermore, we find that compressive strain particularly favors 
a certain rebonded structure over all the others, which explains the atomically smooth and 
periodic structure of the Ge/Si(105) surface 0IH1- 

To show how to obtain structural models for Si(105) in a systematic manner, we start 
from the bulk-truncated structure and attempt to reconstruct the surface in such a way that 
each atom has at most one unsatisfied (dangling) bond after reconstruction. We now present 
in detail the structures of Si(105) with single- and double- height steps. The bulk-terminated 
Si(105), given in Fig. ^a), consists of Si(001) terraces of width 5a/4 separated by steps of 
monatomic height a/4, where a is the lattice constant of Si (a = 5.43A). Terraces can be 
reconstructed by forming short rows of dimers; because the dimer rows are oriented at 45° 
angles with respect to the direction of the steps, every other atom on the step edges must be 
eliminated in order to lower the number of dangling bonds (db). Depending on the relative 
position of the step-edge atoms that are eliminated, there are two distinct ways to achieve 
surface reconstruction. With the notation adopted in Fig. ^b), one possibility is to eliminate 
atoms "1" on edges A and atoms "2" on edges B, and to then form dimers on the surface, 
as indicated by the dotted lines in Fig. Wi c )- This is the model originally proposed by Mo et 
al. ^01) and later named PD (paired dimers) [5J- Because the atoms on the terraces do not 
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rebond at the step edges, we call this model SU (single- height, unrebonded) [3]. 

Another way to achieve single-height reconstruction is by eliminating atoms "1" on 
edges A and atoms "3" on edges B (Fig. ^b)) and then by creating bonds between the 
remaining surface atoms as indicated by the solid and dotted line segments in Fig. ^d). 
This model, which we call SR (single-height rebonded) [Hj, was proposed by Khor and Das 
Sarma and has recently been shown to appear on the side facets of the Ge quantum dots 
0ElEj- We note that there are two different types of bonds on this surface: the usual dimer 
bonds (dotted lines), and the bridging bonds (solid lines). The bridging bonds join a two- 
coordinated atom and a three coordinated one, leaving the former with only one dangling bond 
and fully saturating the latter. Because the bridging bonds of the SR model are stretched, this 
reconstruction is strongly stabilized by the compressive strains present in Ge films deposited 
on Si surfaces [SJ E] - Such rebonded models can also be found for structures that have 
double-height steps, as discussed below. 

We now focus on the unreconstructed (105) surface with double-height steps shown in 
Fig. |2fa)-(b). Since the step height has doubled, the width of the terraces must also be 
doubled, in order to preserve the overall surface orientation. We found several models of 
double-height structures of Si(105), with different terrace structures (e.g. (\/2 x 1) or (2 x 1) 
reconstruction) and atomic bonding at the steps (e.g. rebonded or unrebonded), as explained 
below. 

The simplest model based on (\/2 x 1) terraces is illustrated in Fig. [21(c), where no atoms 
are eliminated and dimers are formed as indicated by the dotted lines. Since no rebonding 
is present, we call this the DU (double-height unrebonded) structure. If we allow for the 
rebonding of atoms on the lower terrace, then all the atoms on the step edges (denoted by C 
in Fig.[2^b)) must be removed, so that any surface atom would have fewer than two dangling 
bonds after reconstruction. The remaining atoms are then bonded as indicated by the solid 
and dotted lines in Fig. Efd). Like the SR model, there are two types of surface bonds, the 
dimers and the bridges; we name this the DR (double-height rebonded) structure. 

In the case of (2 x 1) reconstructed terraces, the dimer rows are oriented at 45° angles 
with the step edge. In analogy to the single-height case, we can eliminate every other atom on 
the step edges (rows C in Fig. Of a)). The elimination of atoms on consecutive terraces can be 
done in-phase or out-of-phase, which leads to structures with different periodic lengths in the 
[501] direction, as shown in Fig. EKb) and (c). Since both of these models involve rebonding 
at the step edges, we call them the DR1 and DR2 structures. From the DR1 structure we 
can obtain the unrebonded model of Zhao et al. j^j by removing another atom from each unit 
cell, as illustrated in Fig.|3Kd); we label this the DU1 structure, to distinguish it from the DU 
model in Fig. |2J 

We have computed the surface energy of the structural models shown in Figs. 11131 using 
two empirical models for atomic interactions, namely the Stillinger- Weber ^2] an d the Tersoff 
[T3~] potentials. While the empirical potentials provide a reasonable description of stepped 
Si(001) surfaces (refer, for example, to ^lj), they are not able to capture the tilting (buckling) 
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of the dimers at the surface, which constitutes an important way of surface relaxation for 
Si(105). Further, it is precisely the tilting of the dimers that determines the major features of 
the STM images and helps in the identification of the atomic structure of the surface El • 
In order to capture the dimer tilting, we have used the charge self-consistent tight-binding 
method of Wang et al. |15| , which accurately predicts the energy ordering of several dimer- 
tilted Si(OOl) structures |TJj]. With this method, for each of the Si(105) structural models 
described above, the total energy of the atoms in the simulation cell exhibits many local 
minima, and we search for the lowest energy structures by using a combination of molecular 
dynamics simulations and annealing. 

Our results are summarized in Tabled where the surface energies of the Si(105) models 
computed using empirical ^3 E] and tight-binding [E] potentials are given. The table also 
contains the number of dangling bonds per unit area for each model. We note that both 
the Stillinger- Weber (SW) and the Tersoff (T3) potentials yield large energy penalties for 
the dangling bonds on the surface, giving an energy ordering similar to that predicted by 
bond counting. On the other hand, the tight-binding (TB) description of atomic interactions 
allows the structures with large numbers of dbs per area (SU, DU and DU1) to relax via 
dimer buckling, leading to an entirely different ordering of the reconstructions. It can be 
seen from Tabled that the surface energies of the (105) models are spread over an interval of 
8-12 meV/A 2 when the empirical potentials are employed. In the case of tight-binding, this 
interval is only ~4meV/A 2 , due to a stronger relaxation of the unrebonded structures. 

For all the potentials that have been used, we find that the SR model has the lowest 
surface energy among all the reconstructions considered; this is in agreement with previous 
work [HI El El where only the SU and SR models are compared. A closer look at the TB 
values in Table Q shows that the SR, DU and DR2 models have energies that fall within 
1 me V/A 2 of the energy of SU, indicating a near-degeneracy of the lowest energy surface. 
Furthermore, our molecular dynamics simulations show that there are other local minima 
(with different bond bucklings) in the same energy interval. Due to this high degeneracy, we 
propose that several models can be present simultaneously on the Si(105) surface, which can 
explain the key features observed in experiments - disorder and roughness. Because energy 
differences are small, the entropy associated with the spatial distribution of different (105) 
unit cells would be important even at low temperatures, explaining the lack of 2D-periodicity 
of the STM images El; the atomic-scale roughness observed in the STM images can be 
generated by a random arrangement of the single- and double-height structures that are close 
in energy. The proposal can be tested experimentally by imaging (zooming in) different areas 
or "patches" of the large-scale Si(105) samples. Preliminary work along these lines has been 
reported in 0. 

We have also examined the strain dependence of the different surface reconstructions. 
By calculating the energies of all the model structures for three values of an applied equibiaxial 
strain (-1%, 0% +1%), we find that the near-degeneracy of the Si(105) surface can be removed 
when the surface is subjected to a compressive state of strain. As illustrated in Fig. even 
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a compressive strain which is as small as -1% further stabilizes the SR model over the other 
models. This result is supported by the recent observations of Fujikawa et al, who showed 
that the initially rough, unstrained, Si(105) surface becomes smooth after being subjected to 
a compressive mismatch strain through the deposition of 3 monolayers of Ge 

The origin of the strain-dependence of the surface energy lies in the arrangement and 
type of the atomic bonds around the step edges: if a structure contains significant rebond- 
ing (i.e. has bridging bonds that are stretched compared to the bulk bonds), compressive 
stresses tend to lower its surface energy. Indeed, this trend emerges from the data presented 
in Fig. all rebonded structures show a decrease of their energy (relative to SU) in compres- 
sion. In contrast, the energy gap between the unrebonded structures (DU and SU, DU1 and 
SU) remains almost constant for the range of strains investigated here. This finding is fully 
consistent with our previous work on the formation energies of (100) steps |3J, as well as with 
the results of Refs. El El on Ge/Si(105). 

In summary, we have constructed a set of structural models for Si(105) and analyzed 
their stability using empirical |12[ I13j and tight-binding J^J potentials. Our study shows 
that the presence of single- and double-height reconstructions on Si(105) can explain the 
experimentally observed |5] atomic roughness and disorder of this surface. Three double- 
stepped reconstructions (DU, DR1, DR2) were found to have lower surface energies than 
the double-step model proposed in jU]. The atomic bonding at the step edge determines the 
strain-dependence of the surface reconstructions and leads to the strain-induced stabilization 
of the single-height rebonded (SR) structure. Future experiments on strained Si(105) surfaces 
(produced, for example, by bending) would be invaluable in gaining further insight into the 
evolution of surface roughness as a function of strain. The models presented here may also 
serve as building blocks for other structures, for example, the quenched (Ix4)-Si(105) observed 
in[3J. 
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Figure 1: Reconstructions of Si(105) with single-height steps and (001) terraces, (a) Side view 
and (b) top view of the bulk-truncated (105) surface. Atoms are colored according to their 
number of dangling bonds (db) before reconstruction: red=2d&, green=lcf&, and blue=0<ife. 
The (105) unit cells are marked by rectangles in figs, (b)-(d). The single-height unrebonded 
(SU) and single- height rebonded (SR) models are shown in (c) and (d); these structures are 
obtained by eliminating atoms "1" and "2", and atoms "1" and "3", respectively (refer to fig. 
(b)). The atoms that are removed are shown as open circles in figs, (c)-(d). The remaining 
atoms are bonded as indicated by black dotted lines (dimer bonds) and purple solid lines 
(bridging bonds). 
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Figure 2: Reconstructions of Si(105) with double-height steps and (y/2 x l)-Si(OOl) terraces, 
(a) Side view and (b) top view of the bulk truncated (105) surface with double-height steps. 
Atoms are colored according to the number of dangling bonds before reconstruction as ex- 
plained in Fig. ^ The (105) unit cells are marked by rectangles in figs, (b)-(d). The double- 
height unrebonded (DU) and the rebonded (DR) structures are shown in fig. (c) and (d), 
respectively. The DR structure is obtained after elimination of all the atoms at the step-edges 
(rows C in (b), open circles in (d)), followed by dimerization (black dotted lines) and rebond- 
ing (purple solid lines) . The thick arrows in figs, (c)-(d) represent the unit vectors of the 
(\/2 x l)-reconstructed terraces. 
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Figure 3: Reconstructions of Si(105) with double-height steps and (2 x l)-Si(OOl) terraces. 
Fig. (a) shows the top view of the bulk-truncated (105) surface with double-height steps. 
Atoms are colored according to the number of dangling bonds before reconstruction as ex- 
plained in Fig. ^ Since the dimer rows on the terraces are oriented at 45° with respect to 
the step edges C, every other atom on the step edges is eliminated upon reconstruction. The 
elimination can proceed in-phase (atoms "1" on all step edges) or out-of-phase (atoms "1" 
on a given step edge and atoms "2" on the next edge). Since rebonding at the step edges is 
present in both models (purple solid lines), we label them DR1 (b) and DR2 (c). The atoms 
that have been removed to achieve the reconstruction are shown as open circles in figs, (c)-(d). 
The (105) unit cells are marked by rectangles with the dimensions 2a x a\/6.5 for DR1 and 
2a x 2a v 7 6.5 for DR2. The unrebonded model DU1 presented in Ref. |Hj can be obtained by 
removing one more atom from the DR1 unit cell as shown in fig. (d). 
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Figure 4: Surface energies of different Si(105) structures measured with respect to the surface 
energy of the SU model, computed with the tight-binding method JH] f° r three values of 
applied equibiaxial strain: +1% (orange), 0%(yellow) and -l%(green). While there are three 
models (SR, DU, DR2) that have surface energies within ~lmeV/A 2 of the surface energy of 
SU at zero strain, a small amount of compressive strain (-1%) removes this near-degeneracy 
and strongly stabilizes the SR model over all the other models. 
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Table 1: Surface energy of Si(105) reconstructions calculated using the Stillinger- Weber po- 
tential (SW) ^21) the Tersoff potential (T3) ^3], and the self-consistent tight-binding method 
of Wang et al. (TB) |15| . The last column indicates the number of dangling bonds (db) per 
surface area expressed in units of a 2 \/6.5, where a is the bulk lattice constant of Si. 
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